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The ATLAS experiment has recently presented its search results for final states 
containing jets and/or b-jet(s) and missing transverse momentum, corresponding to 
^ \ an integrated luminosity of 165 pb~^. We employ this data to constrain a class of 



^ 



super symmetric SU{A)c x SU{2)l x SU{2)r models with t — b — r Yukawa unification, 

in which the gluino is the next to lightest supersymmetric particle (NLSP). The 
cn I NLSP gluino is slightly (~10-30%) heavier than the the LSP dark matter neutralino, 

and it primarily decays into the latter and a quark-antiquark pair or gluon. We find 
(— I I that NLSP gluino masses below ~ 300 GeV are excluded by the ATLAS data. For 

Mh| LSP neutralino mass ~ 200 — 300 GeV and fi > 0, where fj, is the coefficient of the 

C^l MSSM Higgs bilinear term, the LHC constraints in some cases on the spin-dependent 

r~| I (spin-independent) neutralino-nucleon cross section are significantly more stringent 

than the expected bounds from IceCube DeepCore (Xenon IT/SuperCDMS). For 
^.^ , A* < 0, this also holds for the spin-dependent cross sections. 
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I. Introduction 

Low scale supersyminetry, augmented by an unbroken Z^ matter (R-) parity, largely 
overcomes the gauge hierarchy problem encountered in the Standard Model (SM) and also 
provides a compelling cold dark matter candidate. In the mSUGRA/constrained minimal 
supersymmetric model (CMSSM) [l|] , as well as in many other realistic models, the lightest 
neutralino (LSP) is stable J2| with a relic density that is compatible with the WMAP dark 
matter measurements [3|. However, the small annihilation cross section of a pure bino LSP 
with mass of around 100^ GeV does not permit one to easily reproduce the required relic 
dark matter abundance 



An interesting scenario which enhances the bino annihilation cross section is bino-gluino 
co-annihilation. In this case the bino and the relevant NLSP gluino (where NLSP stands for 
next to lightest supersymmetric particle) are sufficiently close together in mass, such that 
the ensuing co-annihilation processes in the early universe allow one to reproduce the desired 
bino relic density. This scenario is not possible in the CMSSM, but it has been implemented 
in models with non-universal gaugino masses [5[, and in a class of (third family) Yukawa 
unified models p-lSJ. The collider sig natures of the gluino co-annihilation scenario have 
recently been discussed in Refs. [ol, llO|. 

The ATLAS and CMS experiments at -y/s = 7 TeV LHC have previously presented their 
search results for low-energy supersymmetry corresponding to an inte grat ed luminosity of 



35 pb~ pjj, 112|, which was recently updated by ATLAS to 165 pb~ pj[. The successful 
launch of the LHC and a flurry of supersymmetry related papers from the ATLAS and 
CMS collaborations provides a strong impetus to explore regions of the MSSM parameter 
space not covered by the minimal version (CMSSM/mSUGRA). In this paper we study the 
constraints and implications of recent LHC data on some well- motivated NLSP gluino models 
induced by gaugino mass non- universality and t — b — T Yukawa coupling unification imposed 
at Mqjjt- The underlying symmetry group we consider is SU{A)c x SU{2)l x SU{2)ji |14l |. 
With the NLSP gluino and LSP neutralino having nearly degenerate masses, the chargino 
as well as leptons are absent in the gluino cascade decay. Also, the jets and missing energy 
from NLSP gluino decay are much softer due to the small mass difference between the NLSP 
and LSP. Thus, the conventional search strategy with same-sign chargino signature does not 
work here, and the usual requirement of large Pt jet and missing transverse momentum 
makes the event selection harder to implement. The LHC constraints on the NLSP gluino 
mass turn out to be significantly less restrictive than the recent 1 TeV or so mass bound on 
the gluino mass which, among other things, assume an essentially 'massless' neutralino. 

The paper is organized as follows. In section II we briefly summarize the NLSP gluino 
scenario with t — b — t Yukawa unification and neutralino (essentially bino) dark matter. 
We also discuss the NLSP gluino decay modes and outline the selection cuts employed by 
the ATLAS collaboration. The results of two classes of NLSP gluino models constrained 
by the LHC data are presented together with a few benchmark points in section HI. Our 
conclusions are summarized in section IV. 
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TABLE L Sunmiary of selection cuts and 95% C.L. upper limits on effective cross section for non- 
SM processes for signal region SI, S2, S3 with 165 pb~^ luminosity, and region b with 35 pb~^ 
luminosity, following ATLAS data analyses [Ij, [iSj] . 



II. NLSP Gluino and ATLAS Selection Cuts 



As mentioned earlier, the gluino-bino co-annihilation scenario requires the gluino to be 
NLSP in the sparticle spectrum, and to be nearly degenerate in mass with the bino LSP. 
The mass difference between the two should be |5| 



\, ""' < 20%. 



(1) 



In the framework of minimal supergravity, this feature clearly requires non-universal gaugino 
masses at Mqut- In particular, a partial unified model given by SU{A)f. x SU{2)i x SU{2)r 
(4-2-2) group structure provides solutions to this scenario. Non-universal asymptotic gaugino 
masses are naturally accommodated in the supersymmetric 4-2-2 model and have recently 
been investigated in Refs. |6|-l8|. With the SM hypercharge in 4-2-2 given by y = y2/5(i? — 
L) + y3/5/3i?, one has the asymptotic relation between the three gaugino masses 

Ml = \m2 + ^Mg, (2) 



where Mi, M2 and M3 denote the asymptotic gaugino masses of t/(l)y, SU{2)i x SU{2)r 
and S'f/(3)c respectively. Assuming that charged fermions of the third family acquires mass 
solely from a single (1,2,2) representation in 4-2-2 leads to the Yukawa unification condition 

at Mgut [li 



yt = yb = Vt 



VDi 



(3) 



It has been shown that t — b — r Yukawa unification can yield relatively light gluinos (< 1 
TeV) 0,0. 

In order to implement radiative electroweak breaking consistent with Yukawa unification, 
the soft mass terms of the two Higgs doublets must be non-universal at Mqut, such that 



the fundamental parameters in this class of models are 

mo, mH^,mHa, M2, M3, Aq, tan/3, sign(/i). (4) 

Here ttiq is the universal soft mass of sfermions, Aq is the universal trilinear scalar coupling, 
tan (3 is the ratio of the vacuum expectation values (VEVs) of the two MSSM Higgs doublets, 
and /i is the MSSM bilinear Higgs mass parameter. The software package ISAJET 7.80 [l7| 
was employed in Refs. J6|-l8| to scan over the relevant parameter space, including renormal- 
ization group evolution of gauge and Yukawa couplings and all soft parameters, as well as the 
computation of the physical masses of all particles. A large number of relevant phenomeno- 
logical constraints such as BR{Bs -^ ju+/i") 0, BR{b -^ s-f) [igj, BR{Bu ^ ru) ^, 



A{g — 2)^ [20I], WMAP relic density J3|], LEP II bound on the lightest Higgs and all the 



sparticle mass bounds 2l| are also implemented. The degree Yukawa of unification is quan- 



tified by the parameter R [22, 6, 7] 

max(2/t, yb, Vr) ,.. 

R = —— . (5) 

We shall require that i? < 1, so that t — b — r Yukawa unification holds at 10% level or 
better. Note that the NLSP gluino scenario with nearly-degenerate gluino and bino masses 
can be realized in 4-2-2 models for both /i > and /x < p-l8|. 

Because of the mass degenerate feature in Eq. ([T]), the NLSP gluino essentially decays 
into colored SM particles such as the gluon octet or a quark-antiquark pair, and the color 
singlet LSP Xv 

9 ^ qqXi,bbx1,gXi, (6) 

where q(q) denotes the first two generation quark (antiquark). The three-body decay g — )■ 
qqXi, bbxi proceeds through an off-shell squark exchange, while the two-body decay g — )■ gxi 
involves a loop diagram containing squarks and quarks. The partial widths of these two decay 
channels are given by |23l . 
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f{x) = 1 + 2x - 8x^ + 18x^ - 18x^ + 8x^ - 2x^ - x' 






- Uxnnx^ + 12x^(1 + x^)lnx^ (9) 

Here Nib, Niw and NiHy, respectively denote the bino, wino and Higgsino components 
of the LSP neutralino Xv Generally, the three-body decays will be suppressed if the scalar 
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FIG. 1. cj X acceptance vs. Mg with horizontal line as the 95% C.L. upper limits on effective 
cross section for non-SM processes for signal region SI (top left), S2 (top right), S3 (bottom left), 
b (bottom right) in the 4-2-2 framework with ^ < 0. Blue regions correspond to models with 
Yukawa unification {R < 1.1). NLSP gluino models form a subset of these and are represented by 
green points. 



masses are too large, or by phase space if the mass difference between g and Xi (AM = Mg — 
M^o) is too small. Assuming M^o ~ (9(250), together with the co-annihilation requirement 
in Eq. ([1]), one has the mass difference AM ^ 50 GeV. Also, for large tan/3, a large 
bottom Yukawa yt naturally leaves the bottom squark (sbottom) to be the lightest squark, of 
C(TeV). With AM ~ 50 GeV and O(TeV) sbottom, g -^ bbxi decay often dominates. One 
can see this feature from Fig. 1 in Ref. [10|, which shows the dependence of the gluino decay 
branching fraction in the Mg — M^ plane for the 4-2-2 model with /i < 0. The NLSP gluino 
decay is therefore sensitive to signals with multi-jets plus missing energy, and relatively more 
to final states with &-jets. 

The ATLAS and CMS collaborations have previously reported data in terms of events 
containing large missing transverse momentum and jets (with or without b-jets) in ^/s = 7 
TeV proton-proton collisions, corresponding to an integrated luminosity of 35 pb~ ^. No 



excess above the Standard Model (SM) background expectation was observed ll|, ll2]- More 



recently, the ATLAS experiment has considered multi-jets plus missing energy events, with 
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FIG. 2. Mg vs. M^o for models with Yukawa unification, NLSP gluino and tliose excluded models 
by ATLAS region SI, S2, S3, b in the 4-2-2 framework with /x < 0. Blue regions correspond to 
models with Yukawa unification (R < 1.1). Green regions are a subset with NLSP gluino, and 
models excluded by LHC data are in red color. 
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TABLE IL Number of excluded 4-2-2 models with Yukawa unification {R < 1.1) and NLSP gluino 
for fi < 0. Also shown are the exclusion fraction by individual signal regions SI, S2, S3, b, and by 
combined SI, S2, S3, b. 



an integrated luminosity of 165 pb~^ ISj. With more strict selection cuts, new lower bounds 
on non-SM cross-sections that are 30 times more stringent than from the 2010 data have been 
obtained. This analyses can also be employed, as we show here, to find useful constraints 
on NLSP gluino models with nearly degenerate gluino and LSP neutralino masses. 

Note that gluino masses below 725 GeV are excluded at the 95% confidence level in 
simplified models containing only squarks of the first two generations, gluino and "massless" 



LSP neutralino 13|. In this case, with the gluino and squarks much heavier than LSP 
neutralino, the large mass difference results in highly energetic jets and large missing energy. 
With nearly-degenerate NLSP gluino and LSP neutralino, however, the jets from gluino 
decay and missing energy are softer, and fewer events with the same gluino mass would pass 
the same selection cuts. We therefore expect that the upper limit on the excluded gluino 
mass for degenerate NLSP gluino and LSP neutralino scenarios would be correspondingly 
lower. 

The CMS analysis has stated less stringent constraints than ATLAS for low-energy su- 



persymmetry search [25|, |26|, and so we utilize the requirements used by ATLAS in our 
studies below. In the updated analysis for multi-jets and missing energy, the events are 
classified into 3 regions "SI", "S2" and "S3", where SI, S2, S3 requires at least 2, 3, 4 



jets respectively [13|. The second class of analysis requires at least one heavy fiavor 6-jet 



candidate in final states 121], denoted by "b" in the following. To simulate similar data, we 



generate all supersymmetric 2 — )■ 2 events and include parton showering and hadronization 



using Pythia [27[, and then forward them to fast detector simulation PGS-4 [28| to simulate 
the important detector effects. The 6-tagging efficiency and mis-tagging rate in PGS-4 are 
based on the Technical Design Reports of ATLAS, and we use the default values in our 
analysis. We further follow the same ATLAS selection cuts for SI, S2, S3 and b. The cut 
requirements are summarized in Tabled where A(j){'p^^^^,ji^2,3) is the smallest azimuthal 
separation between the 0t direction and the three leading jets, and wT-g// is the scalar sum 
oi ^T and the transverse momenta of the highest px jets (up to two for region SI, three 
for region S2 and four for regions S3 and b respectively). The 95% CL. upper limits on 
effective cross section (cross-section times acceptance) for non-Standard Model (SM) pro- 
cesses for signal region SI, S2, S3, b are also showed in Table [B Following Ref. |26| we apply 
a X acceptance > cXexp as exclusion requirement for each model, where a is the relevant total 
cross-section and the acceptance is the ratio of signal events after and before selection cuts 
which reflects the effects of experimental efficiency. 



III. LHC Constraints on NLSP Gluino and Neutralino Dark Matter 

A. t — b — T Yukawa Unification with /i < 

In Refs. p, lZI , the supersymmetric 4-2-2 models with t — b — r Yukawa unification are 
studied for positive p] and negative [TJ] signs of the MSSM parameter fi. The SU{2)l 
gaugino mass M2 was chosen to have the same sign as /i in order to remain consistent with 
the [g — 2)n measurement. This is because the supersymmetric contribution to {g — 2)^ 
is proportional to fiM2. In this section we first consider the ATLAS constraints on 4-2-2 
models with /i < 0. In this case, the finite threshold correction to the Yukawa coupling yi, 



involving the gluino has the desired negative sign. Namely 29 



.„susY-finitc _ _93_ /iMgtan/? _yf_ Mttan/3 
^^ ^ 127r2 M2 ^ 327r2 M^ ' ^ ^ 

where g^ is the strong gauge coupling. At is the stop trilinear coupling, and Mb ~ (M^ -|- 
M^J/2, Mt ^ (Mj2 + /i)/2. For the desired Yukawa unification (^ 10% or better), one 
obtains a wide range of acceptable gluino masses, namely C(200) ^ Mg ^ C(IOOO) GeV as 
shown in Fig. 5 of Ref. [7|. In particular, for relatively light gluinos, Yukawa unification 
is compatible with the gluino-bino co-annihilation mechanism and requires near-degenerate 
NLSP gluino and LSP neutralino masses. 

To study the LHC constraints on this class of models, we generate about half a million 
models by scanning the parameter space |7| and finally obtain 5420 models after applying 
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FIG. 3. asi (left panel) and asn (right panel) vs. M^o in 4-2-2 models with Yukawa unification, 
NLSP gluino, and ^ < 0. The excluded region is denoted in red. The current limits from CDMS-II, 
XENONIOO, SuperK and IceCube and future projected sensitivities from XENONIT, SuperCDMS 
and IceCube DeepCore are also shown. 
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FIG. 4. Mg vs. M^o for models with Yukawa unification (R < 1.1, in blue), NLSP gluino (subset, 
in green) and those excluded by ATLAS region SI, S2, S3, b (in red), for fi > 0. 



the various experimental constraints listed in section HIl Out of these, 3945 models have ac- 
ceptable Yukawa unification [R < 1.1), and in 3807 of these models gluino is the NLSP. The 
region in which the NLSP gluino and LSP neutralino are nearly mass degenerate corresponds 
to 250 GeV < Mg < 300 GeV. 

In Fig. [T] we show ax acceptance vs. Mg for 4-2-2 models with Yukawa unification, 
using the ATLAS regions SI, S2, S3 and b. The near-degenerate NLSP-LSP points are 
also specified and actually overlap with the Yukawa unified points in the low gluino mass 
region. One can see that near-degenerate NLSP-LSP models with Mg < 300 GeV are 



a. 



10- 



-10 



10-" 
lO-ii^ 



H > 0-Region SI, S2, S3, b 




CDMS II 

XtNON Joir 



"••;:;.,XENON 11 
SupiTCDMS~"~"-"--"-"--- 



R < 1.1 
NLSP li 

Excluded 



10 20 30 50 100 

M^; (GeV) 






J3 
.a 



10° 

10-2 

10--' 

10-'' 

io-« 

10- 



Q 

Vi 

b 



10- 
10- 
10- 
10- 



10-2" 

200 300 400 10 



H > 0-Region SI, S2, S3, b 



R < 1.1 
NLSP SI 

Excluded 



Super-K 



20 30 50 100 

M^; (GeV) 




200 300 400 



FIG. 5. a SI (left panel) and gsd (right panel) vs. M^o. Color scheme is same as in Fig. [3l The 
current limits from CDMS-II, XENONIOO, SuperK and IceCube and future anticipated bounds 
from XENONIT, SuperCDMS and IceCube DeepCore are also shown. 



essentially excluded. To display this perhaps more clearly, in Table HT] we outline the number 
of excluded NLSP models with Yukawa unification and the excluded fraction for these models 
by individual signal regions SI, S2, S3, b and combined SI, S2, S3, b. Among the three 
regions SI, S2, S3 of multi-jets plus missing energy final states, region S3 is the weakest 
for constraining NLSP gluino because it requires four jets with px > 40 GeV. However, the 
softest jet from a pair of NLSP gluinos more likely cannot have pt more than about 20 GeV. 
Therefore, a greater number of events do not pass the selection cuts of region S3 compared 
with region SI and S2. Furthermore, as we expect, region b with 6-jets in the final states 
excludes a significant fraction of NLSP gluino models, although we employ the early LHC 
data with 35 pb~^ integrated luminosity. This is because the decay g -^ bbxi is dominant 
in most of the NLSP gluino region, which makes the NLSP gluino models more sensitive to 
multi-6 jets signature. 

In Fig. [2] [Mg — M^o plane), we display (in red color) the excluded models which have 
Yukawa unification and NLSP gluino. One can see that heavier gluinos with Mg ^ 500 GeV 
are consistent with Yukawa unification, but being fairly massive, they survive the current 
LHC constraint. 

It is important to see the implications of LHC data on direct and indirect dark matter 
detection in this class of Yukawa unified models with NLSP gluino. In Fig. [3] we display this 
by plotting the spin-independent and spin-dependent WIMP-nucleon scattering cross-section 
asi (left panel) and asn (right panel) vs. M^o. A significant region around M^o ^ 200 GeV 
is excluded by LHC data, although it is allowed by CDMS-II, XENONIOO, SuperK and 
IceCube experiments. This excluded region will be tested in the future by XENON IT 
and SuperCDMS, but the region lies about three orders of magnitude below the expected 
IceCube DeepCore bound. 
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TABLE III. LHC limits on three NLSP gluino benchmark points that satisfy ah the experimental 
constraints described in Section |TI1 Various selection cuts from the four signal regions, namely, 
SI, S2, S3 and b exclude point 1, whereas point 2 is excluded by the first two regions. Point 3 is 
allowed by all four signal regions. 
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B. t — b — T Yukawa Unification with fi > 

With /i > 0, the gluino contribution to 5y^°'*'^ is positive, so that the contribution from 
the chargino loop must be negative and sufficiently large in order to overcome this. In this 
scenario lower gluino masses and larger values of At and M^ in Eq. ( !T0|) are favored. All 
realistic NLSP gluino models compatible with the WMAP dark matter constraint give rise in 
this case to gluino masses in the range 220 GeV ^ Mg ^ 400 GeV, with the LSP neutralino 
closely degenerate in mass. Also, because of the large At and Mf, values, in this scenario one 
of the stops is usually the lightest squark, with the sbottom relatively heavier than in /i < 
case. Thus, the three-body decay g — )■ bbxl through an off-shell sbottom is suppressed, so 
that the constraint from fe-jets in the final states is less stringent. We start with about 1 
million models and obtain 17942 models which survive the low-energy experiments. Out of 
these, about 400 models display acceptable Yukawa unification {R < 1.1) and contain NLSP 
gluino. Note that the constraint from (g — 2)^ is ignored in generating these models [7|. 

After applying the ATLAS selection cuts listed in Table H] a significant number of models 
are excluded. We show this in Fig.Hlin the Mg — M^o plane. The subset of models with NLSP 
gluino which overlaps with Yukawa unification is also specified. NLSP gluino masses below 
about 250 — 300 GeV are essentially excluded. In Fig. [5] we display the spin-independent 
and spin-dependent WIMP-nucleon scattering cross-section asi (left panel) and aso (right 
panel) vs. M^o. One can see a significant region around M^o ^ 250 GeV is excluded 
by the LHC data although it is allowed by CDMS-II, XENONIOO, SuperK and IceCube 
experiments. Indeed, some parts of the excluded parameter space lie beyond the reach of 
future experiments such as XENON IT, SuperCDMS and IceCube DeepCore. 

Finally, in Table UTTl we present three characteristic benchmark points with NLSP gluino, 
dark matter neutralino and very acceptable t-b-r Yukawa unification. Points 1 and 2, 
with gluino masses close to 300 GeV are excluded by the selection cuts listed in the table. 
However, point 3 with NLSP gluino mass close to 450 GeV is compatible with the data. 



IV. Summary 

Inspired by the recent LHC search of final states containing jets and/or b-jet and miss- 
ing transverse momentum, corresponding to an integrated luminosity of 165 pb~^, we have 
explored its ramifications for supersymmetric SU{4:)c x SU(2)l x SU{2)ii models which dis- 
play t-b-r Yukawa unification at 10% level or better, contain NLSP gluino, and possess LSP 
neutralino dark matter. The NLSP gluino primarily decays into the LSP neutralino and a 
gluon or quark-antiquark pair, thus allowing us to exploit this LHC data. For /i < we 
generate about 4000 models, from an initial half a million models, which satisfy the above 
criteria of Yukawa unification, NLSP gluino, and neutralino dark matter, after imposing con- 
straints from all previous experiments (except LHC). The corresponding number of models 
for /i > is around 400. We next show that for closely mass degenerate NLSP gluino and 
LSP neutralino, models with NLSP glunio masses below 300 GeV or so are largely excluded 
by the LHC data. The LHC implications for spin-dependent and spin-independent LSP 
neutralino-nucleon cross sections are also explored. Regions of the parameter space, some 
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lying well below the much anticipated future bounds from IceCube DeepCore and Xenon 
IT and Super CDMS, are already excluded by utilizing the LHC data. 
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